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Abstract. Context. Using Particle-In-Cell simulations i.e. in the kinetic plasma description, the discovery of a 
new mechanism of parallel electric field generation was recently reported. 

Aims. We show that the electric field generation parallel to the uniform unperturbed magnetic field can be 
obtained in a much simpler framework using the ideal magnetohydrodynamics (MHD) description. 
Methods. We solve numerically ideal, 2.5D, MHD equations in Cartesian coordinates, with a plasma beta of 0.0001 
starting from the equilibrium that mimics a footpoint of a large curvature radius solar coronal loop or a polar 
region plume. On top of such an equilibrium, a purely Alfvenic, linearly polarised, plane wave is launched. 
Results. In ideal MHD the electric field parallel to the uniform unperturbed magnetic field appears due to fast 
magnetosonic waves which are generated by the interaction of weakly non-linear Alfven waves with the transverse 
density inhomogeneity. In the context of the coronal heating problem a new two stage mechanism of plasma 
heating is presented by putting emphasis, first, on the generation of parallel electric fields within an ideal MHD 
description directly, rather than focusing on the enhanced dissipation mechanisms of the Alfven waves and, second, 
dissipation of these parallel electric fields via kinetic effects. It is shown that a single Alfven wave harmonic with 
frequency v — 7 Hz and longitudinal wavelength Aa = 0.63 Mm, for a putative Alfven speed of 4328 km s~^, 
the generated parallel electric field could account for 10% of the necessary coronal heating requirement. It is also 
shown that the amplitude of the generated parallel electric field exceeds the Dreicer electric field by about four 
orders of magnitude, which implies realisation of the run-away regime with associated electron acceleration. 
Conclusions. We conjecture that wide spectrum (10~^ — 10^ Hz) Alfven waves, based on the observationally 
constrained spectrum, could provide the necessary coronal heating requirement. The exact amount of energy 
that could be deposited by such waves through our mechanism of parallel electric field generation can only be 
calculated once a more complete parametric study is done. Thus, the "theoretical spectrum" of the energy stored 
in parallel electric fields versus frequency needs to be obtained. 
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1. Introduction 

The coronal heating problem, the puzzle of what main- 
tains the solar corona 200 times hotter than the 

photosphere, is one of the main outstanding ques- 

2005) for 
significant 

amount of work has been done in the context of heat- 
ing of open magnetic structures in the solar corona 
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tions in solar physics (see e.g. [Tsiklauri 
a recent brief review on the subject). A 
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Historically, all phase mixing studies have been per- 
formed in the Magnetohydrodynamic (MHD) approxima- 



tion; however, since the transverse scales in the Alfven 
wave collapse progressively to z ero, the MHD approxirna. - 
tion is inevitably violated. Thus. lTsiklauri et al.l l)2005allJ ) 
studied the phase mixing effect in the kinetic regime, using 
Particle-In-Cell simulations, i.e. beyond an MHD approx- 
imation, where a new mechanism for the acceleration of 
electrons due to the generation of a parallel electric field in 
the solar coronal context was discovered. This mechanism 
has important implications for various space and labora- 
tory plasmas, e.g. the coronal heating problem and accel- 
eration of the solar wind. It turns out that in the magne- 
tosphere a similar parallel electric field generation mech- 
anism in transversely inhomogen eous plas mas has been 
report ed ijCxenot et al.ll2004l Il999l) . See also lMottez et alJ 
and references therein. 
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This new mechanism occurs when an Alfven wave 
moves along the field in a plasma with a transverse density 
inhomogeneity. The progressive distortion of the Alfven 
wave front due to differences of local Alfven speed then 
generates the parallel electric field. In this work we show 
that the electric field generation parallel to the uniform 
unperturbed magnetic field can be obtained in a much 
simpler framework using an ideal MHD description, i.e. 
without resorting to complicated wave particle interac- 
tion effects such as ion polarisation drift and resulting 
space charge separation, which seems to be the fundamen- 
tal cause of electron acceleration. See also the Discussion 
section for a more detailed account of what led us to con- 
sider the MHD approximation for parallel electric field 
generation. 

In this paper, we also explore the implications of this 
effect for the coronal heating prob lem. Prelimina ry re- 
sults have been reported elsewhere l|Tsiklaurill2006() . The 
present study explores the importance of MHD versus 
kinetic effects in solar plasmas, which has recently at- 
tracted considerable attention due to apparent difficul- 
ties in the coronal heating problem, as well as the natural 
tendency of wave heating models to proceed to progres- 
sively small spatial scales. Also, recent observations have 
indi cated further sub-structuring of the coronal morphol- 
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ogy l)McEwaii fc DeMoorte]l200fi() . which is relevant to our 
earlier work l|Tsiklaurill2005|) . 

2. The model, rationale and main results 

Unlike previous studies I Tsikla^iri et al.l l2fl05albl . 
iGenot et all \2004. Il999() . here we use an ideal MHD 
description of the problem. We solve numerically ideal, 
2.5D, MHD equations in Cartesian coordinates, with 
a plasma beta of 0.0001 starting from the following 
equilibrium configuration: A uniform magnetic field Bq 
in the z— direction penetrates plasma with the density 
inhomogeneity across the x— direction, which varies 
according to 

p{x) = po [1 + 2 (tanh(x + 10) + tanh(-.T + 10))] . (1) 

This means that the plasma density increases from some 
reference background value of po, which in our case was 
fixed at po = 2 X lO^/irrip g cm~^ (with a molecular weight 
of /i = 1.27 c orresponding to the solar coronal conditions 



Hi'^He^lOil. lAschwandenl l|2004) and Trip being the pro- 
ton mass), to 5yOo- Such a density profile across the mag- 
netic field has steep gradients with a half-width of 3 Mm 
around x ~ ±10 Mm and is essentially flat elsewhere. 
Such a structure mimics e.g. the footpoint of a large cur- 
vature radius solar coronal loop or a polar region plume 
with the ratio of the density inhomogeneity scale and the 
loop/plume radius of 0.3, wh ich is the median va lue of 
the observed range 0.15 - 0.5 (|Goossens et al.ll2002l) . The 
above quoted values of coordinate x are in dimensionless 
units. We use the following usual normalisation Bx,y.z = 
BoBx,y,z, {x,y,z) = ao{x,y,z), t = {a(i/c\)t. Here we fix 
Bq to 100 G, and hence the dimensional Alfven speed. 
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Fig. 1. Dimensionless density (Eq.(l)), solid line, and 
Alfven speed, dashed line, profiles across the uniform un- 
perturbed magnetic field (i.e. along x-coordinate) which 
is used as an equilibrium configuration in our model of a 
footpoint of a solar coronal loop or a polar region plume. 



c\ = Ba/^JWp^ turns out to be 4328 km s"^=0.0144 c 
(c is the speed of light). The reference length ao was fixed 
to 1 Mm, i.e. the dimensionless time of unity corresponds 
to 0.2311 s. We usually omit the bar on the top of nor- 
malised quantities, hence when numbers are quoted with- 
out units, we refer to dimensionless units as defined above. 
The dimensionless Alfven speed (normalised to cj^) is then 
ca{x) = 1/^1 + 2 (tanh(j; + 10) + tanh(-a: + 10)). The 
simulation domain spans from — 40Mm to 40 Mm in both 
X— and z— directions with the density ramp as described 
above mimicking a footpoint fragment of a solar coronal 
loop or a polar region plume. Our initial equilibrium is 
depicted in Fig. 1. 

The initial conditions for the numerical simulation are 
By — Acos{kz) and Vy = —CA{x)By at t = 0, which 
means that a purely Alfvenic, linearly polarised, plane 
wave is launched travelling in the direction of positive 
zs. The rest of the physical quantities, Vx and Bx (which 
would be components of fast magnetosonic waves if the 
medium were totally homogeneous) and Vz and Bz (the 
analogs of slow magnetosonic waves) are initially set to 
zero. The plasma temperature is varied as the inverse of 
Eq.(l) so that the total pressure always remains constant. 
Boundary conditions used in our simulations are periodic 
along the z- and the zero gradient along the x-coordinates. 
We fixed the amplitude of the Alfven wave A at 0.05 
throughout. This choice makes the Alfven wave weakly 
non-linear. Our motivation for this value was t wo-fold : 
First, this was the value used bv lTsiklauri et al.1 l 2005b() . 
iGenot et al.l ( 2Q^ , and as one of our goals here is to show 
that the parallel electric field gene ration is also possible in 
the MHD hmit, we use this value. iTsiklauri et alJ ll2005bl) 
only used one value oi A ~ 0.05. while iGenot et alJ 
|200i) used several values of A and showed that weak 
non-linearity (at A ~ 0.05), in addition to the trans- 
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verse density inhomogeneity, of course, is a key factor 
that facilitates parallel electric field generation. Smaller 
As reduced the effect of parallel electric field generation. 
Second, the observed values of the Alfven waves at heights 
oi R = I.OARq ~ 28 Mm are about 50 km s~^ , which for a 
typical Alfven speed of 100 km s~^ makes A equal to 0.05. 
Such observations ( s ee e.g . iMoraiJ l)200l|) . iBaneriee et al.1 
l|l998(l . lDovle et all 



1)1998(1 ) indicate that off limb, heavy 
element (e.g. Si VIII) line broadening is beyond the ther- 
mal one, and hence mostly undamped, outward propagat- 
ing Alfven waves travel through the stratified plasma with 
increasing amplitude (due to the stratification). 

All previous AC-type models of coronal heating fo- 
cused on the mechanisms (e.g. phase mixing or resonant 
absorption) that could enhance damping of the Alfven 
waves. However, for the coronal value of shear (Braginskii) 
viscosity, by which Alfven waves damp of about rj — 1 
m^ s~-^, typical dissipation lengths (e-fold decrease of 
Alfven wave amplitude over those lengths) are ~ 1000 
Mm. Invoking the somewhat ad hoc concept of enhanced 
(anomalous) resistivity can reduce the dissipation length 
to the required value of the order of the hydrodynamic 
pressure scale height At ~ 50 Mm. In this light (obser- 
vation of mostly undamped Alfven waves and the inabil- 
ity of classical (Braginskii) viscosity to produce a short 
enough dissipation length), it seems reasonable to focus 
rather on the generation of parallel electric fields which 
would guarantee plasma heating, should the energy den- 
sity of the parallel electric fields be large enough. In both 
cases (AC-type MHD models and our model with paral- 
lel electric field generation) the energy comes from Alfven 
waves. Therefore observing undamped Alfven waves could 
mean two things: (i) the Alfven waves observed through 
Doppler line broadening are low frequency ones, while high 
frequency Alfven waves produce heating, hence they do 
not contribute to the observed line broadening. Clearly 
one cannot observe waves that have already dissipated; 
(ii) stratification (increase of Alfven wave amplitude with 
height) could be balanced by the damping. 

We emphasise parallel electric fields because in the di- 
rection parallel to the magnetic field electrons (and pro- 
tons) are not constrained. In the direction perpendicular 
to the magnetic field, particles are constrained because of 
the large classical conductivity a — 6x 10^^ s^^ (for T — 2 
MK corona) and inhibited momentum transport across 
the magnetic field. The parameter that controls cross-field 
transport (flow) is cr, and not plasma /3 which merely con- 
trols compressibility of the plasma. This is often confused 
in the literature. 

In the ideal, linear MHD limit there are no paral- 
lel electric fields associated with the Alfven wave. In 
the single fluid MHD the equation for the electric field 
can be obtained by differentiating Ohm's law over time 
and then using Maxwell's equations, expressing rjj by E 
llKrall fc TrivelDiecelll97i : 



V^E- 



Ana dE Ana d V x B 



However, in the coronal plasma conductivity is so large 
that a more simple relation (obtained from Eq.(2)) can be 
used: 



E 



V X B 



(3) 



It is clear from the latter equation that the parallel electric 
field is 



VxBy — VyBx 



(4) 



This means that in the considered system can only 
be generated if the initial Alfven wave (Vy,By) is able 
to generate a fast magnetosonic wave (V^jB^), which, in 
turn, is generated due to weak non-linearity and transverse 
de nsity inhomogeneity. 

iNakariakov et al. 1 lll997l) have investigated this pos- 
sibility of growth of fast magnetosonic waves in a sim- 
ilar physical system, in a context other than parallel 
electric field generation. They used a mostly analytical 
approach and focused on the early stages of the sys- 
tem's evolution. Later, the long-term evolution of the 
fast magnetosonic wave generati on was studied nu mer- 
ically in the case of harmonic l|Botha et al.l l2000(l and 
Gaussian (Tsiklau ri et al.ll200l|) Alfvenic initial pertur- 
bations. When fast magnetoacoustic perturbations are 
initially absent and the initial amplitude of the plane 
Alfven wave is small, the subsequent evolution of the wave, 
due to the difference in local Alfven speed across the x- 
coordinate, leads to the distortion of the wave front. This 
leads to the appearance of transverse (with respect to 
the applie d magnetic field) gradien ts, which grow linearly 
with time. INakariakov et al have shown that with 

fairly good accuracy (which is substantiated by our nu- 
merical calculations presented below) the dynamics of the 
fast magnetoacoustic waves can be described by 



c\{x) 



dt 



\ dx"^ dz 



K = 



(5) 



dt 



dt 



(2) 



Moreover, for small Alfven wave amplitudes the back- 
reaction of the generated fast magnetosonic wave on the 
Alfven wave can be neglected and Eq.(5) can be considered 
as a linear wave-like equation with a driver term on the 
right hand side where By is given by the travelling wave 
expression By = A cos (fc(z — CA{x)ty). The main negative 
outcome of the stu dies with the harmon ic l|Botha et alJ 
|200C|) and Gaussian l|Tsiklauri et alJl200l|) Alfvenic initial 
perturbations was that despite the power-law growth in 
time of the driver term, which implies progressive growth 
of fast magnetosonic wave amplitude (at least until Vx 
and Bx reach the same amplitudes as Vy and By when 
the the back-reaction can no longer be neglected, ren- 
dering Eq.(5) invalid), the Vx and Bx amplitudes after 
rapid initial growth tend to satura te due to the d estruc- 
tive wave interference effect ((Tsiklauri et al.lE00l|) . As a 
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self-consistency test, in the next sub-section we present re- 
sults of the numerical simulation when the wavenumber of 
the initial Alfven wave is fc = 1. In dimensional units this 
corresponds to an Alfven wave with frequency (v = 0.7 
Hz), i.e. longitudinal wave-numbers — 6.3 Mm. 

2.1. The case of an Alfven wave with v = 0.7 Hz 

In Fig. 2 we show two snapshots of Vx and Ez (the lat- 
ter was reconstructed using Eq.(4)) for the case of fc = 1, 
i.e with v = 0.7 Hz. The fast magnetosonic wave (Vx) 
and parallel electric field (Ez) are both generated in the 
vicinity of the density gradients x ~ ±10, eventually fill- 
ing the entire density ramp. This means that the gen- 
erated parallel electric fields are confined by the density 
gradients, i.e. the solar coronal loop which the considered 
system tries to mimic after about 20 Alfven time scales 
becomes filled with temporally oscillating parallel electric 
fields. Note that this can be a source of polarisation drift 
of ions if kinetic effects are considered. 

In Fig. (3) we plot the amplitudes of Vx = Vx and 
Ez = Ez which we define as the maxima of absolute val- 
ues of the wave amplitudes along the x ~ ±10 line (which 
track the generated wave amplitudes in the strongest den- 
sity gradient regions) at a given time step. Solid lines 
represent solutions using the linear McCormack code, 
which solves Eq.(5) with the initial conditions as de- 
scribed above, while dash-dotted lines with open sym- 
bols are the solu tions using the non-linear MHD code 
Lare2d ijArber et al. 2001,). We refer to the McCormack 
code as linear because we treat Eq.(5) as a linear wave- 
like equation with a driver term on the right hand side 
where By is given by the travelling wave form By = 
Acos {k{z — CA{x)t)). Here for both runs we use a 2000 x 
2000 spatial grid resolution. Note that data with solid lines 
is plotted with much smaller time steps than those with 
dash-dotted lines and open symbols, which are plotted 
with the much coarser time step of At = 1. This is due to 
the fact that Lare2d uses MPI parallelisation and hence 
tracking time evolution of the solution at a given point 
of the simulation domain (which is split into many parts) 
is a difficult task. The linear McCormack code is serial, 
hence we do not encounter such problems. The two main 
observations are: 

(i) We gather from Fig. (3) that as in the previous cases 
l)Botha et alJl2000llTsiklauri et alJlioOl]) the amplitude of 
Vx initially grows rapidly and then tends to saturate at a 
level of = 0.05^ = 0.0025, as one might expect from 
the weakly non-linear theory (cf. Eq.(5)). 

(ii) The solutions produced by the fully non-linear 
MHD code and the linear McCormack code are simi- 
lar but not identical. This can be explained by the fact 
that the amplitude oi A = 0.05 is now large enough for 
the non-linearity effects to be noticeable. Note that when 
the amplitude was much smaller, A = 0.001, the devi- 
ation s between the two we re much smaller (cf. Fig. (6) 
from iTsiklauri et"al] l|2n0lh ). In the case oi A = 0.05 
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Fig. 3. Time evolution of the amplitudes oiVx = Vx and 
Ez = Ez- Solid lines represent solutions using the linear 
McCormack code, while dash-dotted lines with open sym- 
bols are the solutions using the non-linear code Lare2d. 
Here A = 0.05, fc = 1, = 0.7 Hz, Aa = 6.3 Mm, and 
both runs have 2000 x 2000 spatial grid resolution. 



with Lare2d, weak non-linearity of the phase-mixed Alfven 
wave even causes noticeable deviations in the initial back- 
ground density, which is an effect previously discussed in 
the literature. 

Figs. (2) and ( 3) corroborate the previous results of 
iBotha et all J2£ 



3) c orroborate the previou s 
and ITsiklauri et alJ l|200ll) . 



2.2. The case of an Alfven wave with u = 7 Hz 

Here we present results for the case of large wave-numbers, 
fc = 10, which in dimensional units corresponds to an 
Alfven wave with ly = 7 Hz, and Aa = 0.63 Mm. This 
is a regime not investigated before. Fig. (4) is similar to 
Fig. (2) but now fc = 10, with shaded surface plots given 
instead of intensity plots. This is due to the fact that 
spiky data do not appear clearly using intensity plots. 
We gather frorn this gr aph t hat similarly to the results of 
ITsiklauri etldl l|2005bl) and lCenot et all l|2004l) . the gen- 
erated parallel electric field is quite spiky, but more im- 
portantly, large wavenumbers i.e. short wavelengths now 
are able to significantly increase the amplitudes of both 
the fast magnetosonic waves (Vx) and the parallel elec- 
tric field Ez. This amplitude growth is beyond a simple 
A^ scaling (discussed in the previous sub-section). The 
amplitude growth is presented quantitatively in Fig. (5). 
The amplitude of Vx now attains values of 0.01, unlike for 
moderate fcs. This boost in amplitude growth can be ex- 
plained qualitatively by analysing Eq.(5). The driver term 
(right hand side of Eq.(5)) contains spatial derivatives. 
Thus, large wavenumbers (i.e. stronger spatial gradients) 
seem to boost the values of the driver term which in turn 
yields larger values for the level of saturation of the Vx 
amplitude. In the considered case, Ez now attains values 
of 0.001. 
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Fig. 2. Intensity plots of Vx and Ez at t = 2 and 20 for the case of fc = 1, = 0.7 Hz, Xa = 6.3 Mm. 
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Fig. 5. Time evolution of the amplitudes of Vx = and 
Ez = E'^. Solid lines with stars represent solutions using 
the Lare2d code with 4000 x 4000 resolution, while dash- 
dotted lines with open symbols are the same but with 
2000 X 2000 resolution. Here A = 0.05, yfc = 10, = 7 Hz, 
Xa = 0.63 Mm. 



Since the amplitude of Vx attains a sizable fraction of 
the Alfven wave amplitude rendering weakly non-linear 
theory inapplicable, we do not plot solutions obtained 
from the the linear McCormack code. Instead, to verify 
the convergence of the solution, we plot the results of the 
numerical run with doubled (4000 x 4000) spatial resolu- 
tion. The match seems satisfactory, which validates the 
obtained results. 



2.3. Application of tine results to the coronal heating 
problem 

It has been known for decades ( Kunerus et al.lll98ll) that 
the coronal energy losses that need to be compensated 
by some additional energy input, to keep the solar corona 
at the observed temperatures, are (in units of erg cm^^ 
s^^): 3 X 10^ for the quiet Sun, 8 x 10^ for a coronal hole 
and 10^ for an active region, ^schwandcn ( 2004) makes 
similar estimates for the heating flux per unit area (i.e. in 



erg cm ^ s ^): 

F„ = EhXt = 5 X 10^ 



(lO^cm) ( 



IMK J 



(6) 



10 ® erg cm s ^. 



This yields an estimate 



2 X 10'' erg cm s ^ in an active region with a 



where Eh 

of Fh « 2 X 10^ erg cm^^ s^^ 
typical loop base electron number density of rig = 2 x 10^ 
cm-3 and T = 1 MK. 

The energy density associated with the parallel electric 
field Ez is 



^-^z r 

Ee^-z—, erg cm 



(7) 



where e is the dielectric permitivity of plasma. The latter 
can be conveniently deduced from the express ion for the 
polarisation current llKrall fc TrivelDiecjlll973l) . appendix 

I) 



pc^ dE 

and Maxwell's equation written as 
VxB.-(J + Jp).-J-,-^-. 



(8) 



(9) 
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Fig. 4. Snapshots of Vx and aX t = 2 and 20 for the case oi k = IQ, v = 7 Hz, Xa = 0.63 Mm. 

the fast magnetosonic waves {Vx and Bx) that propagate 
across the magnetic field and associated parallel electric 
fields {Ez ) are generated in density gradients by the Alfven 
waves. Hence, the flux is carried with the Alfven speed) to 
obtain 

Fe = Eec\ = 1.62 X 10^ [ergcm-^s-i] , (14) 

which is « 10 % of the coronal heating requirement esti- 
mate for the same parameters made above using Eq.(6). 
Note that the latter estimate is for a single harmonic with 
frequency 7 Hz (see the Discussion section for details when 
a wide spectrum of Alfven waves is considered). 

3. Dissipation of the generated parallel electric 
fields 

We now discuss how the energy stored in the generated 
parallel electric field is dissipated. We examine the paral- 
lel electric field behaviour at a given point in space as a 
function of time. In Fig. (6) we plot the time evolution of 
Ez at a point [x — 10.68, z = 0) for the case of fc = 1, 
V = 0.7 Hz, and \a = 6.3 Mm. Choice of this x-value is 
such that it captures parallel electric field dynamics at the 
strongest density gradient point (across x). 

We gather from Fig. (6) that Ez is a periodic (sign- 
changing) function that is a mixture of two harmonics 
with frequencies uj — CAk and w = 2c Ak. This is due 
to the fact that E^ is calculated using Eq.(4) where Vy 
and By at fixed spatial points vary in time with frequen- 
cies UJ = CAk, while t he generated Vx and B x vary with 
frequencies to — 2c Ak llNakariakov et al.ll997j) . Because of 
the ideal MHD approximation used in this paper, the gen- 
erated electric field cannot accelerate plasma particles or 



The result is 

s-'-^- (10) 

The latter formula is differen t from the usu al expression 
for the dielectric permitivity ijKrall fc Trivelpiecei.1973,) : 

47rpc^ ,,,, 
£ = l + -g^- (11) 

because the displacement current has been neglected in 
the above treatment, which is a usual assumption made 
in the MHD limit. For the coronal conditions {p — 2 x 
lO^fimp g cm^-^, /i — 1.27, B = 100 Gauss) the second 
term in Eq.(ll) is 4.8048 x 10^ w e > 1 which means that 
for the coronal conditions in the energy density Eq. (7) the 
dominant term is 

l^I^l = \,Vr, (12) 

which is the kinetic e nergy density associated with the 
ExB drift tSturrocJ l)l994j) . page 60). In Fig. (5) we saw 
that electric field amplitude attains a value of « 0.001. 
In order to convert this to dimensional units we use 
= 4328 km s'^ and B = 100 G and Eqs.(3-4) to ob- 
tain Ez ~ (caB/c) X 0.001 = 0.0014 statvolt cm'^ (in 
Gaussian units). Therefore the energy density associated 
with the parallel electric field Ez (From Eq.(7)) is 

Ee^ex 0.0014V(87r) = 3.7471 x lO"'' (13) 

[ergcm^"^] . 

In order to get the heating flux per unit area for a sin- 
gle harmonic with frequency 7 Hz, we multiply the latter 
expression by the Alfven speed of 4328 km s~^ (because 
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Fig. 6. Time evolution of Ez at a point [x = 10.68, z = 0). 
Here fc = 1, i/ = 0.7 Hz, and Xa = 6.3 Mm. 



cause Ohmic heating unless kinetic effects are invoked. Let 
us look at Ohm's law for ideal MHD (Eq.(3)) in more de- 
tail denoting physical quantities unperturbed by the waves 
with subscript and ones associated with the waves by 
a prime; initial equilibrium implies Eq = Vq — with 
Bo 7^ 0. For the perturbed state (with Alfven waves {Vy 
and By) launched which generate fast magnetosonic waves 
(Vx and Bx) ) we have 



E' = -V X {Bo + B')/c. 



(15) 



Note that the projection of E' on the full magnetic field 
(unperturbed Bq plus the waves B') is zero by the defini- 
tion of the cross and scalar product : E' ■ {Bo + B')/\{Bo + 
B')\ = 0. Physically this means that in ideal MHD the 
electric field cannot do any work as it is always perpen- 
dicular to the full (background plus wave) magnetic field. 
However the projection of E' on the unperturbed mag- 
netic field Bo is clearly non-zero 



E' 



Bo 

\Bn 



= Ez = - 



V X B' Bo 



\Bo\ 



(16) 



which exactly coincides with Eq.(4) that was used to cal- 
culate Ez- The crucial next step that is needed to un- 
derstand how the generated electric fields parallel to the 
uniform unperturbed magnetic field dissipate must invoke 
kinetic effects. In our two stage model, in the first stage, 
bulk MHD motion (waves) generate parallel electric fields, 
which as we saw cannot accelerate pa rticles if we d e scribe 
the plasma in the ideal M HD limit. iGenot et all ^^2004^ 
and lTsiklauri et all l|2005b(l showed that when the identi- 
cal system is modelled in the kinetic regime particles are 
accelerated with such parallel fields. Genot et al... 1,2004) 
claimed that electron acceleration is due to the polarisa- 
tion drift. They showed that once the Alfven wave propa- 
gates along the density gradient that is transverse to the 
magnetic field, a parallel electric field is generated due to 
charge separation caused by the polarisation drift asso- 
ciated with the time-varying electric field of the Alfven 



wave. Because the polarisation drift speed is proportional 
to the mass of the particle, it is negligible for electrons, 
hence ions start to drift. This causes charge separation 
(the effect that cannot be treated by a MHD description), 
which results in generation of parallel electric fields, that 
in turn accelerate electrons. As was shown in Fig. (6), our 
parallel electric field is also time- varying, hence at the sec- 
ond (kinetic) stage the electron acceleration can proceed 
in the same manner through ion polarisation drift and 
charge separation. The exact picture of the particle dy- 
namics in this case can no longer be treated with MHD 
and a kinetic description is more r elevant. The frequencie s 
considered in the kinetic regime llTsiklauri et alJl2005bl) . 

= 0.3cjci/(27r) = 4.6 x 10"* Hz and 0.7 - 7 Hz (this pa- 
per) are different, but we clearly saw that the increase 
in frequency results in enhancement of the parallel elec- 
tric field generation. Various effects of wave-particle in- 
teractions will rapidly dampen the parallel electric fields 
on a time scale much shorter than the MHD time scale. 
ICenot et~a l. (2004) clearly demonstrated the role of non- 
linearity and kinetic instabilities in the rapid conversion 
from the initial low frequency electromagnetic regime to a 
high frequency electrostatic one. They identified Buneman 
and weak beam plasma instabilities in their simulations 
(as t hey studied the t i me evo lution of the syst em for longer 
tha,n iTsiklaiiri et abl l|2005hl)V Fig.(ll) from lOenot et "abl 
l|2004|) shows that wave energy is converted into parti- 
cle energy on times scales of lO^uip^ ~ 4 Alfven periods. 
Perhaps no immediate comparison is possible (because of 
the different frequencies involved), but still in our case the 
Alfven period when sufhcient energy is stored in the par- 
allel electric field is 1/(7 Hz) =0.14 s, i.e. wave energy is 
converted into particle energy in a short time. 

Yet another important observat ion can be ma de by es- 
timating the Dreicer electric field l|Dreiceilll959l) . Dreicer 
considered dynamics of electrons under the action of two 
effects: the parallel electric field and friction between elec- 
trons and ions. He noted that the equation describing elec- 
tron dynamics along the magnetic field can be written as 



rUeVd 



eE - Vp/'nieVd 



(17) 



where Vd is the electron drift velocity and Vp^^ is the elec- 
tron collision frequency and dot denotes time derivative. 
When Vd <C Vthermai Eq.(17) in the steady state regime, 
{d/dt — 0) allows us to derive the expression for Spitzer 
resistivity. When Vd > vthermai, the steady state solution 
may not apply. In this case, when the right hand side of 
Eq.(17) is positive i.e. when eE > Vp^^rUeVd, we have elec- 
tron acceleration. This is the so-called run-away regime. 
In simple terms, acceleration due to the parallel electric 
field leads to an increase in Vd] in turn, this leads to a de- 
crease in Vp^^ because Vp^'' oc l/wj^. Thus, when the electric 
field exceeds the critical value Ed — n^e^ In A/(87reQfcBT) 
(the Dreicer electric field, which is quoted here in its SI 
form), faster drift leads to a decrease in electron- ion fric- 
tion, which in turn results in even faster drift and hence 
the run-away regime is reached. Putting coronal values 
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(rie = 2 X 10^ cm-^, T = 1 MK and InA = 17.75) 
in the expression for the Dreicer electric field we obtain 
0.0054 V which in Gaussian units is 1.8 x 10~^ stat- 
volt cm~^. As can be seen from this estimate, the am- 
plitudes of the parallel electric field obtained in this pa- 
per exceeds the Dreicer electric field by about four orders 
of magnitude. This guarantees that the run-away regime 
would take place leading to electron acceleration and fast 
conversion of the generated electric field energy into heat. 
In Eq.(17) E is constant, while our Ez at the point of 
strongest density gradient is time varying. Hence some 
modification of the Dreicer analysis is expected. 

4. Discussion 

After the comment paper bv lMottez et all l|2006(l we came 
to the realisation t hat the electron acceleration seen in 
both series of works l|Tsiklauri et al .l2nn5albl . lc^iMt' et al.1 

is a no n-resonant wave - particle interaction ef- 
fect. In works bv IXsiklanri et all 1 2005albt ) the electron 
thermal speed was vth,e = 0.1c while the Alfven speed in 
the strongest density gradient regions was va = 0.16c; 
this unfortunate coincidence led us to the conclusion 
that the electron acceleration by parallel electric fields 
was affected by the Landau r esonance with the phasc- 
mixed Alfven wave. In works bv lCenot et all l 200^T l9991 
the electron thermal speed was Vth,e = 0.1c while the 
Alfven speed was va = 0.4c because they considered 
a more strongly magnetised plasma applicable to Earth 
magnetospheric conditions. However, the interaction of 
the Alfven wave with a transverse density plasma inho- 
mogeneity when the Landau resonance condition lu = 
fc||CA(a;) is met can be quit e impor tant for the el ectron ac- 
celeration JChaston et al.li2000.. Haseeawa fc Chen,1976) . 
IChaston et alJ l)2000() assert that the electron acceleration 
observed in density cavities in aurorae can be explained 
by the Landau resonan ce of the cold ionospheric electrons 
with the Alfven wave. Illasegawa fc ChenI l)l976|) also es- 
tablished that at resonance the Alfven wave fully con- 
verts into the kinetic Alfven wave with a perpendicular 
wavelength comparable to the ion gyro-radius. We can 
then conjecture that because of kinetic Alfven wave front 
stretching (due to phase mixing, i.e. due to the differences 
in local Alfven speed), this perpendicular component grad- 
ually realigns with the ambient magnetic field and hence 
creates the time varying parallel electric field component. 
This points to the importance of the Landau resonance 
for electron acceleration when the resonance co ndition is 
met. But as seen in works of lGenot et all l|2004. 1999 ), 
even when the resonance condition is not met, electron 
acceleration is still possible. 

There were three main stages that lead to the formu- 
lation of the present model: 

(i) The realisation that the parallel electric field gener- 
ation (and particle acceleration) is a non-resonant wave- 
particle interaction effect lead us to the question: could 
such parallel electric fields be generated in a MHD ap- 
proximation? 



(ii) If one considers non-linear generation of the fast 
magnetosonic waves in the transversely inhomogeneous 
plasma, then E = —{V x B)/c contains a non-zero 
component parallel to the ambient magnetic field Ez = 

-{V,By-VyB,)/c. 

(iii) From prev ious studies l|Botha et al.l l2000l 
iTsiklauri et al.ll200l[l we knew that the fast magnetosonic 
waves {Vx and B^) did not grow to a substantial fraction 
of the Alfven wave amplitude. However, after reproducing 
the old parameter regime (k=l, i.e a frequency of 0.7 Hz), 
the case of k=10, i.e a frequency of 7 Hz was considered, 
which showed that fast magnetosonic waves and in turn a 
parallel electric field were more efficiently generated. 

There are two main issues that need to be discussed 
to address the plausibility of the proposed model: (i) the 
plausibility of parameters used in the model and (ii) the 
relation to the observations. 

First, the parameter space of the problem is quite 
large. The level at which the fast magnetosonic wave {Vx 
and Bx) and hence the parallel electric field {Ez) ampli- 
tudes saturate d epends on several parameters. As indi - 
cated previously ijBotha et al.l2000llTsiklauri et al.l200l|) . 
this level depends on the strength of the transverse den- 
sity gradient and the wavelength (width in the case of a 
Gaussian Alfvenic pulse) of the Alfven wave. In particu- 
lar it was found that the stronger density gradients yield 
lower saturation levels of the fast magnetosonic wave am- 
plitude (due to the fact that destructive wave interference 
starts earlier when the density gradients are strong) , while 
shorter wavelengths (width in the case of the Gaussian 
Alfvenic pulse) of the Alfven wave generate h igher levels 
(e.g. Figs. 11 and 12 in ITsiklauri et alJ l)200l|) ). We have 
not done a full parameter space investigation to demon- 
strate the effect, but instead we fixed the transverse den- 
sity gradient guided by observations. In particular, the 
observed length scale of the density inh omogeneities in 
loops vary between 0.15 and 0.5 loop radii ( Goossens et alJ 
I2OO2I) . In our model the length scale of the density inho- 
mogeneity (half-width) is 3 Mm and the loop radius is 
10 Mm (see Fig. 1) which makes the ratio 0.3. This is 
the median value in the observed range (0.15 ~ 0.5). This 
eliminates one parameter of variability in the parameter 
space. For the Alfven speed, we used a putative value of 
4328 km s^^ and performed two numerical runs for two 
frequencies, 0.7 and 7 Hz. A full investigation should map 
a range of frequencies. Such an analysis is pending, but see 
below for preliminary estimates in the context of coronal 
heating. 

Second, Alfven waves as observed in situ in the so- 
lar wind always appear to be propagating away from the 
Sun and it is therefore natural to assume a solar origin 
for these fluctuations. However, the precise origin in the 
solar atmosphere of the hypothetical source spectrum for 
Alfven waves (turbulence) is unknown, given the impos- 
sibility of remote magnetic field observa tions above the 
chrom osphere-corona transition region (IVelli fc Primeti 
I1997I) . Studies of ion cyclotron resonance heating of the 
solar corona and high speed winds exist which pro- 
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Fig. 7. Plot of the flux carried by Alfven waves from 
10~^ Hz up to a frequency v versus that frequency 
as inferred from t he em pirically constrained spectrum 
of ICranmer et alJ l)l999(l . The solid line corresponds to 
Eq.(18), while the dashed line shows te coronal heating 
requirement for a temperature of 1 MK. 



vide important spectroscopic constr aints on the Alfven 
wave spectrum l|Cranmer et al Although the spec- 

trum can and is observed at distances of 0.3 AU, it 
can be projected back to the base of corona using em- 
pirical constraints, se e e.g. the top line in Fig. 5 from 
ICranmer et all \\m^ (see also t he mo re elaborate model 
of ICranmer h van Balleeoouenl l|2005|) '). Using the latter 
figure we can make the following estimates. Let us look 
at single harmonic, first. At a frequency of 7 Hz (used 
in our simulations), the magnetic energy of Alfvenic fluc- 
tuations is e\]'^'^^ k, lO'' nT^ Hz"-'^. For a single har- 
monic with — 7 Hz this gives for the energy density 
£;(7Hz) ^ i.£;(7Hz)/(8^) ^ 7 ^ ^0-3 GV(87r) « 2.8 x IQ-^ 
erg cm~'^. Surprisingly this semi-observational value is 
quite close to the theoretical value given by Eq.(13). As 
we saw above, such a single harmonic can provide approx- 
imately 10 % of the coronal heating requirement. Next let 
us look at how much energy density is stored in the Alfven 
wave spectr um based on the empir ically guided top line in 
Fig. 5 from ICranmer et al.l ljl999f) . Their spectral energy 
density (which they call "power") is approximated by the 
so-called 1// spectrum, i.e. E,, = 0.6 x 10^ /i^ nT^ Hz"^ 
In proper energy density units this is Ei^ — 2.4 x 10^'^ /v 
erg cm^'^ Hz Thus, the flux carried by Alfven waves 
from say 10^"* Hz up to a frequency v would be 



F 



AW 



E^c^^dv = 1.04 X 10^ ln(iy/10"'^) (18) 



[ergcm^^s^^] . It is instructive to look at this equation 
graphically in Fig. (7): the Alfven wave spectrum from 
1/ = 10^^ Hz up to about a few xlO'^ Hz carries a flux 
that is nearly as much as the coronal heating requirement 
(Eq.(6)) with T — 1 MK. We do not consider higher fre- 
quencies because at about 10"* Hz, ions become resonant 



with circularly polarised Alfven waves and dissipation pro- 
ceeds though Landau resonance - a well-studied mecha- 
nism, but quite different from our non-resonant mecha- 
nism of parallel electric fleld generation. 

There are several possibilities of how this flux carried 
by Alfven waves (fluctuations) is dissipated. If we consider 
the regime of frequencies up to 10"^ Hz, the ion cyclotron 
resonance condition is not met and hence dissipation is 
dominated by the mechanism of parallel electric field dis- 
sipation formulated in this paper. However, at this stage it 
is unclear how much energy could actually be dissipated. 
This is due to the fact that we only have two points, 0.7 
Hz and 7 Hz, in our "theoretical spectrum". As we saw, 
a single Alfven wave harmonic with frequency 7 Hz can 
dissipate enough heat to account for 10% of the coronal 
heating requirement. Eq.(18) shows how much flux is car- 
ried by the Alfven waves, while in order to calculate how 
much of it is actually dissipated depends on what level Fe 
from Eq.(14) will attain for each harmonic. Hence, the 
flux dissipation through our mechanism would be 



E^c\D{v)dv 



(19) 



where D{v) is the "theoretical spectrum" of the energy 
stored in parallel electric flelds. For a 7 Hz single harmonic 
it can be obtained from 



2.4 X l{r^c\D{7Y{i) 



(7 Hz) = Eec\ 



(20) 



with Ee = 3.7471 x lO"'^ [ergcm^^j fj.^^^ Eq.(13) ren- 
dering Z?(7Hz) — 1.56. For a 0.7 Hz single harmonic Ee 
would be different because as can be seen from Fig. (3) 
E^ K i X 10-5 (as opposed to E^ 0.001 for v ^ 7 
Hz). Since Ee scales as El, then D{{).7m) = 1.56 x 
(3 X 10-5/0.001)2 « 1.4 X 10-3. More numerical runs are 
needed to fully map D{v). This is deferred until further 
work is done. 

A clear distinction should be made between our model 
and the ones that use ion cyclotron resonance damping of 
Alfven waves. Low frequency (0.001-0.1 Hz) Alfven waves 
dominate the power spectrum of ffuctuations in the solar 
wind at and beyond 0.3 AU. These waves are able to trans - 
port a signiflcant amount of energy l|Cranmer 
On the contrary, high frequency (10-10000 Hz) Alfven 
waves are known to damp more easily than their low fre- 
quency counterparts, but they are not expected to con- 
tain much power. This is problematic for ion cyclotron 
resonance damping models because Landau resonance of 
ions occurs at high frequencies (few 10** Hz). Our model 
on the contrary is of a non-resonant nature. Thus, it is 
hoped that it can provide enough heating once the fre- 
quency range lO^"' — 10'^ Hz at the base of corona is 
mapped numerically (i.e. once the level attained by the 
parallel electric field amplitudes for each frequency is de- 
termined). The proposed wide spectrum idea for Alfven 
waves is not as " theor etical-rather-th an-demonstrable" as 
the one proposed by iTsiklauri fc Nakariakov, (2iIQi) for 
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the slow magnetosonic waves. One could question the va- 
lidity of the wide spectrum idea for the slow magnetosonic 
waves in coronal loops as predominantly single harmonics 
(with periods of 3 and 5 minutes, etc.) are observed. At 
the same time it is not possible to observe the high fre- 
quency waves that are dissipated already. However, in the 
case of Alfven waves, the presence of a wide spectrum 
with a frequency range of 10~* — 10^ Hz at the base of 
corona, which is actually observed at 0.3 AU, seems to 
be more likely. The possibility of meeting the full coronal 
heating requirement with the wide spectrum Alfven waves 
via the proposed two stage mechanism of parallel electric 
field generation needs further investigation. 
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